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Cross Reference to Related Applications 

Benefit is claimed, under 35 U.S.C. 1 19(e), to the filing date of provisional patents 
applications serial numbers 60/420,090 and 60/419,834 , entitled "Method for analyzing light 
simultaneously of all wavelengths scattered off 2D or 3D micro- or nano-electronic structures", and 
"Method for significantly reducing the oscillatory nature of scattered polarized radiation to better 
enable the solution of the inverse problem of determining the properties of the scatterer^', both filed 
on 10/21/2002 and both are listing Eytan Barouch, Brookline MA; as inventor, for any inventions 
disclosed in the manner provided by 35 U.S.C. 1 12. These provisional applications are expressly 
incorporated herein by reference. 

Field of the Invention 

The thrust of this invention is a new methodology for very fast and accurate of structure and 
composition of objects at wide range of sizes, scattering a wide band of electro-magnetic radiation. 
This invention includes analytical, algorithmic and practical numerical implementation by analyzmg 
the reflected light off of micro- and nano-electronics features for advanced metrology applications, 
thereby enabling the sizes and features of the structures to be determined on location in real time in a 
manufacturing environment. 

Background of the Invention 

This patent application has to do with a method and algorithm that can be used in conjunction 
with scatterometry apparatus. The following background, quoted fi*om the backgroimd of US patent 
#5,963,329, serves fairly well here also as a starting point on the optical measurement of micro and 
nanoelectronic structures: 

"In microelectronics, accurate measurement of feature profiles (i.e. line width, line height, 
space between lines, and sidewall shape) are very important to optimizing device performance and 
chip yield. Measiurements are needed at many steps in manufacturing to assure that critical 
dimensions, line profiles, and feature depths are under control. Historically, measurements have been 
accomplished with the following technologies: 

Optical imaging using the resolving power of optical microscopes and iniage processing to 
measure small features. However, features smaller than the resolving power of the microscope can 
not be measured, nor can the line profile. 
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; Electron-beam imaging, particularly the scanning elec^ 
improves resolution over the optical microscope. However, like optical imaging, top-down SEM 
imaging does not provide profile information. While imaging a cross-sectioned samples does 
provide profile information, cross sectioning is destructive, costly, and labor intensive. In addition 
the signalprocessing used in SEM imaging introduces uncertainties. Furthermore, electron beams 
can daniage the sample, and this is especially the case when sensitive materials such as photoresist 
are imaged. In addition, electron.beam charging can seriously distort measurement signals, causing 
beam position deviations. To overcome this problem, conductive coatings are conmionly used but 
coatings or their removal can damage the sample or effect measurement accuracy. 

. Scanning force microscopy (SFM) utilizes a very small mechanical probe to sense minute 
atomic forces that act very close to the surface of materials, and SFM can provide high-resolution 
topographical information, including line profiles. Unfortunately, SFM is slow and operator- 
intensive, and the quantitative science of the probe tip is uncertain and unreliable. 

Surface-contact mechanical-probe technologies, such as surface profilometers, have been 
used to profile large structures but do not offer the resolution or sensitivity required by today's 
characterization requirements. Furthermore, contact-probes distort surfaces. In addition, mechanical 
stability requirements prohibit the use of probes small enough to acconmiodate the sub micrometer 
sizes commonly measured. To measure the width of a line the location of each edge of its profile ^ 
must be defined, and an arbitrary, qualitative edge detection model is usually used in each of the 
above meaisuring techniques. While this arbitrary measurement point is typically calibrated to a 
cross-section, manufacturing process changes and normal manufacturing process drift can invalidate 
the edge model and introduce significant measurement error. 

, . Furthermore, none of the above techniques can characterize plaharized, buried sub-micron 
structures that commonly are found in microelectronic or micro-machining manufacturing processes. 
Features may be formed that are later covered with a transparent (or semi-transparent) layer of 
material, which is then planarized. If the feature sizes are below the resolution limits of optical 
techniques and the flat top surface prevents SEM and AFM imaging, none of the above well-known 
techniques will work. 
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Two. optical techniques are available for planarized buried structures. Scatterometry has been 
used for characterizing periodic topographic structures in which a light beam (typically a laser) 
illuminates an area to be characterized and the angular distribution of elastically scattered light is 
used.to measure line widths on photomask and silicon wafer gratings. However, scatterometry 
requires a high dejgree of mechanical stability and accurate measurement of the laser beam 
impingement angle on.the grating, and these mechanical limitatioiis contribute to the uncertainty of 
scatterometry results. 

Confocal laser imaging is another optical technique that permits characterization of 
planarized buried structures. Focus on the sample is progressively adjusted with a monochromatic, 
high numerical aperture (NA) optical imaging system, and those focus adjustments provide the line 
profile. The high NA lens has a very small depth of focus, and therefore, a physical profile can be 
obtained by moving through a series of focus settings up or down the features, and recording those 
height variations. However, the technique is limited by the resolution of the laser beam, or the sizie of 
the laser spot on the sample. Furthermore, the use of a high NA objective lens, providing a large 
angle to the incident beam limits the ability to characterize high aspect ratio features that are 
commonly found in semiconductor manufacturing." 

Patent number 5,963,329 went a long way toward improving the above situation. As quoted 
in the summary of that patent: 

"It is therefore an object of the present invention to provide an optical method of determining 
a line profile of a sub-micron structure. It is a feature of the present invention that the line profile is 
determined by comparing measured diffraction radiation intensity from a periodic structure on a 
substrate with modeled radiation intensity predicted from a theoretical line profile. It is a feature of 
the present invention that line profile parameters of the model are varied until the predicted intensity 
converges to the measured intensity. It is a feature of the present invention that a rigorous and 
detailed two-dimensional analytical solution is obtained for the line profile. It is an advantage of the 
present invention that the line profile is reliably, reproducibly, rapidly, and nondestructively, 
obtained. These and other objects, features, and advantages of the invention are accomplished by a 
method of determining the profile of a line comprising the steps of: (a) providing a substrate having 
a repeating structure comprising a plurality of lines, said lines having substantially identical profiles; 
(b) illuminating said repeating structure with radiation wherein said radiation diffracts, said 
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diffracted radiation having an intensity; (c) measuring said int^sity; (d) providing a model structure 
on.a data processing machine,! said model structure comprising a repeating structure on said 
substrate, said model structure comprising a model profile; (e) mathematically predicting a predicted 
diffracted radiation intensity when said model structure is illxmiiriated with said radiation; and (f) 
comparing said predicted intensity with said measured intensity." 

The present invention has to do witii an advanced algorithm that works well with either 
scatterometry or ellipsometry equipment, but in particular shows its highest advantage when used 
with analyzing light gathered in the usual method of scatterometry. 

Now, there is another class of measurement technology, namely, ellipsometry, that we need 
to mention here that does involve the measurement of light reflected off of substrates at an angle, 
and measuring the change in intensity and polarization. By describing some of the operation of this 
measurement procedure, some of the aims/purposes of the present invention will become clearer. In 
a conventional ellipsometry set-up, the angle of incidence and the angle of reflection of the beam, to 
the normal of the substrate, are made as nearly equal to each other as possible. Ellipsometry is well 
known for being able to deduce the complex index of refraction of the substrate, for each wavelength 
of light used in the illumination. More specifically, by shining well coUimated light, at some angle to 
the normal, of a substrate material, and by collecting and measuring the intensity of the reflected 
light at the same mgle to the normal, and by determining the change in polarization, then it is well 
known that one can deduce what the complex index of refraction, n(a))=n(a)) +ik(a)) will be, for 
each wavelength of light for which this measurement is made. If the substrate has a single, vmifonhly 
composed fihn layer on it, and if one knows the substrate's optical properties [c.g.ynfk) ], and if the 
film is transparent, then the real part of the complex index of the film, and the film's thickness, can 
also be deduced. If one has a stack of fihns, then one can also determine the effective complex index 
of refiraction of the stack of fihns, provided other information is supplied. 

In general, in ellipsometry, the angle of incidence to the normal of the substrate, and the 
angle of reflection from the substrate, are made identical For each state of incident radiation (known 
wavelength, known polarization state, and known angle of incidence to the sample), two 
measurements can be made, namely, the intensity of the light scattered into the reflected beam, for 
each of two polarization states that can always be used to characterize the reflected beam. Thus, 
there are two measurements that can be made for each incident ray and polarization state onto the 

5 



sample..Ih.the case ofa flat substrate, without structure and without fihns, then these two 
measurenients can be used to deduce n and k for the substrate. If a fihn exists on the substrate, then n 
and k for the substrate, and n and k for the fihn, and the thickness d, cannot all be deduced, as one 
only has two pieces of information (change in intensity of Ught into the reflected beam for each 
exiting polarization state). Howeyer, by measuring the substrate properties first without the film, 
then supplying the information of n and k of the substrate in the second set of measurements of 
involving the fibn. on the substrate, then, one can reduce the nuihber of unknowns in the coiiiposite 
system. There are still too many unknowns, however, namely, n, ky and d for the film. However, if 
thie film is transparent, so that A=0, or if independently one can measure the absorption properties of 
the fibn, then one can use the two sets of data obtained fi"om the ellipsometry measurement to 
deduce n and d for the fibn. Similar, but more complicated procedures, can be pursued to obtain 
information about stacks of uniform fikns [2]. 

In general ellipsometry practice, there are two other commonly used means for collecting yet 
more information about the planar substrate.and composite fihn structure.The same set of 
measurements just described can be done carried out as a function of wavelength, by simply using a 
broadband light source and filtering the incident light to make it of a narrow wavelength, or quasi- 
monochromatic, character. For each incident light situation of distinct wavelength, one can carry out 
the same measurements just described. Thus, one can measure the change in intensity of each 
polarization state, for more than one wavelength of incident radiation, to thereby gain information 
regarding the complex mdex 8(a)>=w(a>)-/*((D), where «(©) and *(<o) are fimctioiis of the frequency 
CO. Furthermore, by carrying out this same measurement scenario for different angles of incidence for 
the radiation on the substrate plane, then yet another entire multiplicity of measurements can be 
obtained. Thus, for each angle of incidence of Ught, and for each wavelength of light iised in the 
illmnination, there are two measurements that can be obtained in the beam of reflected light, where 
the reflected light angle to the normal, in ellipsometry, is niade equal to the incident angle of the 
directed light. Most of these points will have bearing on the invention to be described here. Thus, 
ellipsometry has to do with measuring the properties of light reflected off at an angle from a 
substrate, with possible layers of vmiform film on the substrate, where the incident and reflected light 
beams have the same angle to the normal to the substrate. Scatterometry has to do with the detection 
of scattered radiation off of a substrate, typically with a periodic structure in or on the substrate, 
where one or more aiigles of scatter from the normal are examined, and where the intensity and 
polarization of the light, for each wavelength, are compared with that of the incident radiation. 
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. : As for the algorithm that is conventionally used in analyzing light obtained via scatterometry 
methods, the typical technique is to follow a rather computationally intensive approach involving 
optimizing the number and sizes of a large number of "slabs"of photoresist to yield a reflectivity 
versus wavelength curve that most closely match the intensity data. Other researchers have also 
made investigations in scatterometry along these lines, or via matching precalculated libraries of 
shapes. The present invention does not have any of these limitations. 

Brief Summary of the Invention 

In this invention we address these considerations in by developing a system of processes 
described in figure 1. The first task is to evaluate the system being examined for general properties, 
type of materials involved and their electromagnetic properties. An important decision is whether 
said.materials contain some inherent non-linear polarization jproperties. The sizes of the features 
(nano, micro or much larger) as well as the range of needed wavelengths (or fi*equencies) to resolve 
and reconstruct the features is evaluated at this initial step. Next are three simultaneous steps: 
material characterization as a function of fi*equency immediately transformed to time domain, a 
"rendering algorithm' characterizing the scattering object with its seed composition and shape and 
spatial computational mesh with its corresponding time steps and stability and precision criteria, so a 
cost function for exit can be defined. The next step consists of inserting the various materials and 
scattering objects into the developed very-fast and accurate forward scattering algorithm using 
hybrid method similar to multi-grid analysis, computing the TE, TM polarization of the reflected . 
waves and then the non-polarized spectrum, to go into the cost function for comparison against the 
target spectrum. Then the structure and composition are updated by various methods (like least 
square or conjugate gradient) to provide an updated composition and shape to be put back into the 
system iteratively. When the cost-tolerance is obtained, the system exits and reports its findings. 

Brief Description of Drawings 

Figure 1 illustrates the invention organizational chart. 

Figure 2, 3 and 4 illustrate the powerful innovation embedded in the module "Findpoles". In 
said figures a comparison between target optical parameters and their computed permittivity function 
is demonstrated, with excellent agreement. 
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Figure 5. illustrates the coinputed arid measured spectra using the module "Fwd". The result 
of "Firidpoles" is applied to a current realistic stack of four layers. 

Figure 6, 7, 8 and 9 display a photolithography state of the art system, comparing 300nm and 
400nm pitches of lines and spaces (L/S). The industrially measured spectra under manufacturing 
conditions are employed as targets for the module "Invscaf and the closest agreement is given. The 
reconstructed feature agrees with the SEM measurements to about .4%. 

Description of the Invention 

This patent application involves a new methodology to solve Material-Maxweirs equations 
for electromagnetic scattered beam of composite lossy substrates and isolated or periodic complex 
and composite objects containing metals as well as absorbing lossy materials. As such, it allows the 
microchip industry to detect and evaluate critical dimension,, affording determination of the 
manufacturability of micro and nanoelectronic striictures. This method of analysis is equally valid 
for all wavelengths of interest and can be done simultaneously in real time on a single processor for 
transverse electric (TE) and for transverse magnetic (TM) polarizations. The Material-Maxweirs 
equations naturally decompose into two separate systems (TE and TM) that can be solved 
simultaneously. Earlier work as described in US patents # 5,963,329, 6,433,878 addresses both facts 
that the range of validity is very Umited and the versatility of shape reconstructidn abiUty- is. limited 
to parallel layers only due to the inefficiency of the RCWP theory. Also, both are limited to one 
wavelength since per calculation since they use an algorithm vaUd only in frequency space. Another 
limiting applicability of the RCWP method is its complete disregard to frequencies in the high and 
low domain, resulting in severe limitations at the deep UV (ultra-violet) range, where most current 
action takes place. 

In this invention we address these considerations in by developing a system of processes 
described in figure 1. The first task is to evaluate the system being examined for general properties, 
type of materials involved and their electromagnetic properties. An important decision is whether 
said materials contain some inherent non-linear polarization properties. The sizes of the features 
(nano, micro or much larger) as well as the range of needed wavelengths (or frequencies) to resolve . 
and reconstruct the features are evaluated at this initial step. Next are three simultaneous steps: 
material characterization as a Amotion of frequency immediately transformed to time domain, a 
'"rendering algorithm' characterizing the scattering object with its seed composition and shape and 
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spatial computational mesh with its corresponding time steps and stability and precision criteria, so. a 
cost: ftm'ction for exit can be defined. The next step consists of inserting the various materials and 
scattering objects into the developed very-fast and accurate forward scattering algorithm using 
hybrid method similar to multi-grid analysis, computing the TE^ TM polarization of the reflected 
waves and then the nonrpolarized spectrum, to go into the. cost function for comparison against the 
target spectrum. Then the structure and composition are updated by various methods (like least 
square or conjugate gradient) to provide an updated composition and shape to be put back into the 
system iteratively. When the cost-tolerance is obtained, the system exits and reports its findings. 

The most critical components ofthe invention are: . 

a. The material composition function is obtained by the use of the module "Firidpoles". This 
module takes a seed table exhibited in Table 2 with the target table exhibited in Table 1 . 
Table 2 expresses the seed as a sum of Debye, I^renz, Lx)renz, Conductivity^^ 
Constant terms. It is then converted to the format jgiven in Table 1 , whose columns are 
the energy in eV, n, k, real permitivity, imaginary permitivity, reflectance and 
wavelength in nm. A complicated least-square code "Findpoles" is turned, on and when 
toleranceis obtained it exits with.the output file displayed in Table 3. There, the inherent 
firequencies of the material are given in terms of 1 /sec; as well as the relaxation 

: . . . . parameters. The strength of each term is dimensidnless and the conductivity is given in; 
. terms ofSiemens. As a by product, *Tindpoles" provides an efficient arid useful ^ 
extraction method of m&ky that satisfies causality restrictions unlike the claim in US . 
patent # 4,905,170, where causahty is violated, rendering it useless for electromagnetic 
calculations in time domain. 

b. The speed of the forward scattering module "Fwd" is largely obtained due to the ability 
. to express the very large domain of substrate material in terins of a very few grid grid- 
cells (-'10). This innovation is achieved by a mathematical formulation depicting the 
large substrate layer as having a "slave" attached lossy layer with equal impedence of the 
substrate material, thus cutting down the most expensive part ofthe computational 
domain substantially, obtaining thousands-fold speed-up. The "slave" layer has poles of 
the same firequency as the substrate material with different coefficients obtained fi*om 
equalizing the impedances. 
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0.0018935 


2.11098 


^.0055021 


0.0340946 


790 



Tablel : Iiiput target for 'Tindpoies' 
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find poles 2. 1 .Ob input data 


1 . 1 . . 


conduct term 


1 


plasma terai 


1 


num lorentz 


5 


num Xlorentz 


1 


nimi debye 


0 


420.00 


84 


791 


nlayers 


4 




min 


initial 


max 


eps_inf 


0.1000000015 


0.8200767636 


7.0000000000 


Lsigma 


0.0001000000 


4.8033914566 


5.0000000000 


omega_pl 


1.0000000000 


6.5119667053 


15.0000000000 


alpha_pl 


0.0000000000 


0.9998378158 


1.0000000000 




S 


0.1000000089 


6.1372642517 


9.3000001907 


omega 


0.4000000060 


1.5478570461 


7.0000000000 


ganmia 


0.0099999998 


0.0100028720 


1.0000000000 




S 


0.0399999991 


2.3346211910 


4.4000000954 


omega 


0.4000000060 


3.5085842609 


13.0000000000 


gamma 


0.0010000000 


0.0379080027 


1.0000000000 




S 


0.0399999991 


0.9923810363 


4.4000000954 


omega 


0.4000000060 


6.6716189384 


12.0000000000 


gamma 


0.0010000000 


0.2266924083 


1.0000000000 




S 


0.0399999991 


1.6832729578 


4.4000000954 


omega 


0.4000000060 


4.3492865562 


9.0000000000 


gamma 


0.0010000000 


0.0379193872 


1.0000000000 




S 


0.0149999997 


4.4830975533 


8.6999998093 


omega 


0.8000000119 


4.0481534004 


8.0000000000 


gamma 


0.0010000000 


0.0641734675 


1.0000000000 




XS 


0.1000000089 


1.8365435600 


2.2999999523 


Xomega 


0.4000000060 


5.2369093895 


12.0000000000 


Xgamma 


0.0099999998 


0.1255621761 


1.0000000000 


Xlambda 


0.0099999998 


0.6012777152 


1.0000000000 



Table 2: Seed input file into "Findpoles' 
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Matfile for input into "Fwd" 


Output of "Findpoles" 


1 eos inf 


0.71569538 


1 siffma 


8 17777292E+12 


1 numlorentz 


5 


1 nuniXlorentz 


0 


1 nuriidebve 

X XX W^XX%*%r T W 


1 ■ ; •■ 






deleps_p 


8.52794349E-01 


omega_p2 


6.39123 152E+31 


delta_p 


6. 1470867 1E+ 14 






deleos d 


1.57363355E+00 


ome&a o2 


2.81803558E+31 


delta, n 


1 51713529E+14 

X X / X*/*^A>^JU^ ' x~ 






del ens n 


1 68789959E+00 


otnesa o2 


9 01482831E+31 


delta n 


2 16371420E+15 






deleps_p 


4.39927959E+00 


omega_p2 


4.27885802E+31 


delta_p 


3.57961 164E+14 






deleps_p 


2.52639604E+00 


omega_p2 


3.41288632E+31 


delta_p 


3.66798193E+14 






delepsD_p 


2.49912262E+01 


tauD_p 


4.95714492E-16 



Table3: Output material file of "Findpoles" as input file into "Fwd" 
The reconstruction of the scattering object is performed with a hybrid method similar to 
the LM algorithm with one fundamental variation. The initial Hessian is obtained by 
localized parameter minimizations, avoiding the extreme inaccuracy involved with 
numerical derivatives. The initially constructed Hessian provides both seed as well as 
vector-directions utilizing a Gram-Schmidt procedure. 
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Figures 2, 3 and 4 illustrate the powerful innovation embedded in the module "Findpoles". In 
said figures a comparison between target optical parameters and their computed permittivity function 
is demonstrated, with excellent agreement. The resuh of "Findpoles" is applied to a current realistic 
stack of four layers and the computed and measured spectra are exhibited in figure 5, using the 
module "Fwd". In figures 6,7,8,9 a photolithography state of the art system is displayed, comparing 
300nm and 400nm pitches of, lines and spaces (L/S). The industrially measured spectra imder 
manufacturing conditions are;employed as targets for the module ""Invscat" and the closest 
agreement is given. The reconstructed feature agrees with the SEM measurements to about .4%. 

The method described in this invention is accomplished by providing a detailed analysis of 
the dielectric function of each material in the structure of interest. The functional form of the 
dielectric function can be expressed in one of two ways. First, it can be given as a sum of a constant 
and an interactive plasma term (non-magnetized) and several poles of Debye, Lorenz double poles 
and Linear-Lorenz double poles. Second, the dielectric function may be expressed in terms of the 
density of states in the corresponding quantum mechanical band, as expressed in an integral equation 
of Cauchy type, i.e., with a singular kernel, that satisfies the fundamental Kramers-Kronig 
relationship and is solved by a combination of analytic and nimierical techniques. An approximate 
analytic solution of the integral equation corresponds to a high level of accuracy to the first method 
just mentioned. The permitivity so described is a representation of the dielectric function in 
fi-equency space, and utilized as such forces solution of the Material-Maxwell equations for each 
wavelength separately, resulting in a very slow execution time as well as severe inaccurate and 
unreliable results. Consequently, despite the patent of ref. 1, IBM has initially rejected this 
technology. Furthermore, the dielectric function in fi'equency space must satisfy several strict 
restrictions in order for it to represent the real physics it attempting to describe. These restrictions 
include the Kramers-Kronig relations, causality, positivity of the index of refiraction n and the 
absorption coefficient k, as well as the property thjat the complex conjugate of the permittivity 
function as a function of firequency co is a function of -oo. The real part of the permittivity function is 
allowed to be positive and negative continuously, presenting another usually ignored difficulty when 
the real part vanishes. When that happens all classical electromagnetic descriptions become invalid 
in fi-equency space and disregarding it for the fi'equency in which it takes place yields absurd results. 
A typical material in which this takes place is Silicon, which is a critical material for the microchip 
industry. The physics behind this occurrence is a change in response of the material &om semi- 
conductor to a temporary plasma containing material. This invention deals in part with overcoming 
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this complication. The time domain dielectric function is the Fourier transform of the permittivity . 
function in frequency space, and as such it is not subjected to the separation of the real and 
imaginary parts. However, there is a high price to pay. The Material-Maxwell equations are no 
longer partial differential equations since they become integral-partial-differehtial (IPD) system, 
which is much harder to handle. An integrial part of this invention is the incorporation of the properly 
characterized material functions into the IPD system allowing a very rapid evaluation of the forward 
scattering for all wavelengths contained in the source, simultaneously. . 

For the microchip industry, the Material-Maxwell solver of this invention uses a few 
nanometer grids. However, the featurie that designates the target needs to be characterized up to a 
fraction of a nanometer. Thus, a rendering algorithm has been developed that divides each grid cell 
into an arbitrary localized refined mesh where the various weights on the original grid cells are 
modified according to the materials included in that cell. For example, if each boundary cell is 
divided into a 10x10 grid, the CD accuracy thus obtained is 0.5 nm. This technique, recently 
developed and unpublished, is employed to compute the inverse problem of determining the 
scattering object shape and material composition from the knowledge of the incoming and the 
scattered beams. It also allows alignment determination of alignment markers and shapes. The speed 
and accuracy of these algorithms would enable in-situ and on-line metrology for realistic systems. . 

This invention is equally applicable to 2D&3D micro and hanoelectronic structures. By the 
development of this method, inspection of OPC (optical proximity corrections) can be obtained in a 
manufacturing environment, while taking into account aUgnment and twist location structures. 
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